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ABSTRACT 

Laser altimeter measurements of forest stands with distinct structures and compositions have been acquired at 532 nm 
(green) and 1064 nm (near- infrared) wavelengths and parallel and perpendicular polarization states using the Slope 
Imaging Multi-polarization Photon Counting Lidar (SIMPL). The micropulse, single photon ranging measurement 
approach employed by SIMPL provides canopy structure measurements with high vertical and spatial resolution. Using 
a height distribution analysis method adapted from conventional, 1064 nm, full-waveform lidar remote sensing, the 
sensitivity of two parameters commonly used for above-ground biomass estimation are compared as a function of 
wavelength. The results for the height of median energy (HOME) and canopy cover are for the most part very similar, 
indicating biomass estimations using lidars operating at green and near-infrared wavelengths will yield comparable 
estimates. The expected detection of increasing depolarization with depth into the canopies due to volume multiple- 
scattering was not observed, possibly due to the small laser footprint and the small detector field of view used in the 
SIMPL instrument. The results of this work provide pathfinder information for NASA’s ICESat-2 mission that will 
employ a 532 nm, micropulse, photon counting laser altimeter. 
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1. INTRODUCTION 

Over the past 15 years lidar remote sensing has proven to be a highly effective method for characterization of forest 
canopy height and vertical structure 11 ' 21 . Numerous papers have shown that this information is a strong predictor of 
above-ground biomass, providing an important means to quantify terrestrial ecosystem carbon storage and through time 
characterize sources and sinks of atmospheric C0 2 . A conventional laser altimeter measurement approach used for this 
purpose is medium-footprint, full-waveform recording typically done in the near-infrared (NIR) P1 . Micropulse, photon 
counting, typically done at green wavelengths, is a more recent approach with improved vertical resolution 131 . In this 
work we assess observations of forest canopy height and structure using this new approach, evaluating the results at 
near-infrared and green wavelengths and two polarization states. Specifically, we evaluate cumulative height 
distributions of detected photons adapting a method developed for conventional waveform measurements. We compare 
the height of median energy (HOME) and canopy cover, standard waveform metrics of canopy structure used in biomass 
estimation, to determine if they are wavelength sensitive and therefore a potential source of difference in biomass results 
obtained by NIR and green measurement approaches. We also assess if wavelength and polarization differences in the 
height distributions can provide information on the scattering of the pulsed laser energy within canopies that could yield 
information on the type and spatial organization of foliage, information potentially useful for differentiating tree species. 
Data used in this study was acquired using the Slope Imaging Multi-polarization Photon-counting Lidar (SIMPL) flown 
over three forested sites having a diversity of canopy vertical structures. 

1.1 Medium-footprint, full-waveform lidar 

Full-waveform systems include three laser altimeters developed by NASA, the airborne Scanning Lidar Imager of 
Canopies by Echo Recovery (SLICER) 141 , the airborne Laser Vegetation Imaging Sensor (LVIS) 151 and the Geoscience 
Laser Altimeter System (GLAS) 16 ' 71 flown aboard the Ice Cloud and Land Elevation Satellite (ICESat) 18 ' 91 . Those systems 
utilize a single pulse analog measurement approach depicted schematically in Figure 1. The systems operate at 1064 nm 
(NIR) with pulse widths in the range 7 to 10 nsec full-width half-max (FWHM), equivalent to 1 to 1.5 m. The pulses 
illuminate footprints on the Earth’s surface with diameters varying between 10 m for SLICER up to 70 m for GLAS. 



Reflected laser energy is detected using a Silicon Avalanche Photodiode (Si:APD) that converts received energy to an 
output voltage time series that is digitized, nominally with 1 nsec (15 cm) sampling, providing a lidar waveform. The 
waveform is a measure of the within-footprint vertical distribution of canopy surfaces (foliage, stems, branches) and the 
underlying ground where illuminated by the laser pulse 141 . High pulse energies are used in order to detect sufficient 
signal (1,000s of photons) above the Si:APD noise floor operating from medium altitude airborne platforms (5 mJ from 
10 km for LVIS) to spaceflight altitudes (50 mJ from 600 km for ICESat). This signal level is necessary to fully 
characterize the canopy structure and to detect the ground surface beneath dense, closed canopies, both of which are 
necessary to determine canopy height. Because of the high pulse energies, the pulse rate must be low for a practical laser 
transmitter implementation (500 Hz for LVIS, 40 Hz for ICESat) thereby limiting their spatial sampling and resolution. 
LVIS acquires contiguous footprints, nominally 20 m in diameter, across a several km wide swath. ICESat samples along 
a single profile, nominally with 50 m diameter footprints spaced apart by 170 m. In addition to limited spatial resolution, 
the relatively broad pulse width acts as a low pass filter, thereby smoothing the waveform and limiting the vertical 
resolution of canopy features. This also broadens the return from the ground and reduces its amplitude thus making its 
detection more difficult 131 . 
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Fig. 1 . Single pulse analog and micropulse photon counting measurement approaches for characterization of forest canopy 
height and vertical structure, adapted from' 31 . 

One method to characterize canopy structure and estimate above-ground biomass using conventional lidar waveforms 
utilizes cumulative distributions of the received energy (Figure 2) f4,ltM21 . Comparison to field plots has shown that a 
multi-parameter, linear regression using several waveform-derived parameters is a strong predictor of field-measured 
above-ground biomass. An early approach 110111 used four parameters: the canopy height (distance from the center of the 
ground return to the highest detected signal above a threshold), the height of median energy (HOME, the distance from 
the ground return to the 50 th %), the ratio of the two heights, and the received energy of the ground return divided by that 
of the canopy (a measure of canopy cover, the percentage of ground covered by a vertical projection of the canopy 
surfaces). Other studies use the ground received energy divided by the total energy as the measure of cover 141 . This is the 
criteria used in this study. HOME is sensitive to changes in both the vertical arrangement of canopy elements and the 
canopy cover 1101 . As observed by a lidar, the cover is a function of both plant area and the retro-reflectance of the ground 
relative to that of the canopy at the laser wavelength 141 , with higher relative ground retro-reflectance resulting in a lower 
cover result. The enhanced-amplitude, “hot spot” retro-reflectance at 0° phase angle is the applicable reflectance for lidar 
instruments because they have parallel illumination and view angles 141 . 

1.2 Micropulse, photon counting lidar 

To address the analog waveform limitations of low detector sensitivity and restricted vertical and spatial resolution an 
alternative measurement approach has been developed. Micropulse, photon counting (Figure 1) uses a high repetition 
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Fig. 2. A lidar waveform and the derived cumulative distribution indicating percentile heights of received energy, including 
the 50 th % height of median energy (HOME) (adapted from https://lvis.gsfc.nasa.gov). 

rate laser transmitter, nominally at 10 kHz with a ~1 nsec (15 cm) pulse width and low pulse energy, ~ 1/1,000 that of a 
traditional system with comparable performance 13 ' 13 ' 161 . Typically a photomultiple tube (PMT) in Geiger mode with high 
sensitivity at visible wavelengths is used to detect the arrival of back-scattered single photons at 532 nm that are 
transmitted by frequency doubling a 1064 nm laser. To increase the count rate of single photons per laser fire a PMT 
array with multiple detector channels is often used. The arrival time of the single photons is nominally resolved to 0.1 
nsec (1.5 cm) using an event timer. The round-trip travel time between the laser fire and photon detection event provides 
the ranging distance to the reflecting target. By accumulating the detected photons for a sequence of laser fires into a 
range histogram a measurement analogous to a wavefonn is generated 131 . Because the micropulse width is — 1 Ox 
natTOwer and the timing precision is ~10x better, a measurement with substantially better vertical resolution than a 
conventional wavefonn is achieved. A key benefit is that returns from canopy layers and the ground are narrower and 
therefore higher in amplitude making their detection more reliable. Fewer detected photons are therefore required to 
achieve a measurement of comparable sensitivity. In photon counting systems the noise floor of the detectors is 
essentially zero and the dominant noise source becomes solar background illumination, which can hinder the detection of 
weak returns. To minimize the noise rate a small detector instantaneous field of view (IFOV) only a few times larger 
than the footprint is used, and narrow band spectral filtering further reduces the noise. NASA’s ICESat-2 mission, a 
follow-on to ICESat planned for launch in 2016, will use this measurement approach 1 ' 71 . Whereas ICESat employed a 
single laser beam, ICESat-2 will measure six profiles in a “push broom” configuration using a diffractive optical element 
(DOE) to split the output of a single laser into six beams. The beams are configured for observing ice sheet elevation 
change by precisely repeating measurement tracks through time. The planned design will have three tracks separated by 
3 km for spatial coverage, with each track comprised of two profiles spaced cross-track by < 100 m in order to separate 
elevation differences due to slope from actual elevation change. Each profile will be sampled with 10 m diameter laser 
footprints and 0.7 m spacing between footprints, providing substantially better spatial resolution than that of ICESat. 


2. SLOPE IMAGING MULTI-POLARIZATION PHOTON-COUNTING LIDAR 


SIMPL is an airborne, micropulse, photon counting laser altimeter developed through the NASA Earth Science 
Technology Office Instrument Incubator Program. It incorporates a variety of advanced technologies in order to 
demonstrate measurement approaches of potential benefit for improved, more efficient spaceflight laser altimeter 
missions including ICESat-2. Instrument details are documented elsewhere 1 ' 81 . In summary, SIMPL incorporates beam 
splitting of a micropulse laser, single-photon ranging and polarimetry technologies at green (532 nm) and NIR (1064 nm) 
wavelengths. The basic measurement approach and the micropulse transmitter characteristics are similar to that 
described above. Like ICESat-2, the measurement configuration is a “push-broom” but with four beams spaced apart 
equally perpendicular to the flight direction. However, because use of a DOE would separate the transmit paths of the 
two colors, a dichroic beam splitter based on bi-refringent calcite crystals was developed to preserve co-aligned colors. 
The splitter also maintains transmit pulses that are plane polarized. For each beam the received energy is separated into 




the two colors and each color is further separated into energy with polarization states parallel (||) and perpendicular (J_) 
to the transmitted polarization, yielding four channels per beam and a total of 16 channels. Traditional PMT detectors 
are not used because they are insensitive atNIR wavelengths. Instead, for each channel SIMPL employs a Single Photon 
Counting Module (SPCM) 1191 that has sensitivity at both visible and NIR wavelengths. The probability of detecting (PD) 
a photon per laser fire is kept below 30% for each channel, by adjusting mechanical irises and flight altitude, in order to 
minimize biasing the ranging distance and depolarization measurements that can be introduced at high detection rates 
due to the long, ~ 50 nsec (7.5 m) dead-time of the SPCM detectors 118,191 . 

SIMPL’s measurements provide information about surface elevation, roughness and slope as well as scattering properties 
used to differentiate surface properties and types. By measuring the depolarization of the plane-polarized transmit pulses 
at the two wavelengths SIMPL has a unique capability to differentiate surface types, building upon prior non-ranging 
measurements of laser depolarization at 355, 532 and 1064 nm [2WI ', By adding the ranging component, SIMPL provides 
a means to investigate the vertical distribution of optical scattering properties to better understand the interaction of 
pulsed laser energy with surface targets. The depolarization ratio (J_ / ||) is sensitive to the proportions of specular 
reflection and surface and volume scattering, and is a function of wavelength. Upon single scattering from a surface the 
reflected energy remains parallel to that of the plane-polarized laser pulse. As multiple scattering increases the fraction 
of the energy converted to the perpendicular state increases as does the depolarization ratio. Prior work using data 
acquired over Lake Erie ice cover demonstrated differentiation of open water, snow and several ice facies as well as 
differences in the depth of laser pulse penetration 1181 , results of relevance for planned ICESat-2 measurements of sea ice 
thickness 1171 . For canopies, the degree of depolarization is a function of wavelength and the reflectance phase function, 
transmission, size and spacing of foliage and woody components and the specular character of waxy broadleaf coatings, 
if present 1201 . Ground depolarization is a function of wavelength, the type of ground-cover and its retro-reflectance and 
roughness. 

2.1 Forest canopy data collection 

SIMPL data was collected over four forest sites using the NASA Lear-25 operated by Glenn Research Center in late July 
and early August, 2010: the Smithsonian Environmental Research Center (SERC), MD where the forest structure has 
been documented by extensive ground observations 1221 ; two locations within the New Jersey Pine Barrens (PB) (Cedar 
Bridge and Silas Little) centered on USDA Forest Inventory and Analysis plots; and a site in the Huron National Forest 
(HNF), MI. Short segments of the data acquired on one flight line at each site are shown in Figures 3-5. The SERC and 
Pine Barrens sites are being used by the ICESat-2 project for measurement validation and algorithm development. At 
each of these sites 15 km data segments were collected on 6 flight lines that repeated lines used earlier in 2010 by a 
swath mapping, 532 nm, micropulse, photon counting lidar operated by Sigma Space Corp. Two lines were flown and 
then repeated at the FINF site that was selected because it consists of single-species, managed pine stands that provided 
distinct changes in forest height and structure. At the nominal flight altitude of 2.3 km above ground level (AGL) at 
SERC, the spacing between the beams was 5 m, the footprint diameter was 0.2 m and the detector IFOV was 0.6 m. 
These dimensions are small because SIMPL’s optical design was predicated on eventual flight on the International Space 
Station at a nominal altitude of 360 km, from which the beam spacing and footprint size would be 800 m and 30 m, 
respectively. The nominal altitudes at the PB and HNF sites were 2.0 km and 1.8 km AGL, respectively. At the laser 
fire rate of 1 1.4 kHz and a nominal aircraft ground speed of 100 m/sec and PD between 5% and 15%, the typical density 
of signal photons detected was between 5 and 15 per m. SIMPL’s micropulse transmitter normally has a 2.7 nsec 
FWHM after pulse occurring 8 nsec (1.2 m) after the primary pulse. The amplitude of the after pulse at 1064 nm is 
about 20% that of the main pulse. The amplitude is reduced to about 5% of the main pulse at 532 nm due to the process 
of frequency doubling. The after pulse was present during the PB and HNF flights. The after pulse amplitude was 
greatly reduced during the SERC flight, inferred to be due to alteration of the laser operating state by high pre-flight 
ambient temperatures in the aircraft. 

2.2 Results 

Plotting ranging distance as a function of acquisition time for each detected photon yields single-photon “point clouds” 
for each of SIMPL’s 16 channels, providing a detailed depiction of the two dimensional structure of the forest stands and 
a profde of the underlying ground topography. Figures 3-5 show examples for the third beam’s 532 nm || channel. The 
point density for the SERC example appears lower than the other two sites because a shorter flight segment is depicted, 
thereby expanding its horizontal axis. Undulations in the range to the ground are a combination of topographic relief and 
changes in aircraft altitude and attitude. Three stands at each site were selected to derive and evaluate range histograms. 
Stands with distinct canopy height, structure and cover were selected based on visual examination of the point cloud 



profiles and associated imagery. At SERC the stand types are: A) —115 year old mixed broad-leaf hardwoods with high 
canopy cover, B) ~70 broad-leaf tulip poplar and sweet gum with moderate canopy cover, and C) 1 5 year loblolly pine, 
selectively logged in 2007, with low canopy cover 1231 . At PB specific information on the selected stands is not available. 
The region has naturally occurring mixed pine and deciduous cover. At HRN the stand types are: A) a recent jack pine 
plantation with high canopy cover, B) a red pine plantation established in 1942 with moderate canopy cover, and C) the 
1942 red pine plantation which has had low thinning applied to reduce under-story tree density 1241 . Distinct changes in 
forest height and vertical structure are especially apparent in the HRN point cloud at boundaries between the stands of 
different age and silviculture management practices (Figure 5). To the right of stand A is a very young jack pine 
plantation stand into which is cut trapezoidal gaps as a part of a program of habitat maintenance for the endangered 
Kirtland's Warbler. The gaps in the short stature pines are observed in the point cloud. To the left of stand B is a recent 
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Fig. 3. Results for SERC, MD data acquired on 07/26/10. Top: SIMPL profile track (red) displayed on an aerial 

photograph from Google Earth acquired on 10/01/08. The location along the track for the three selected stands is 
indicated by the black segments, each of which are ~50 m in length. At this scale the small cross-track separation 
between the four beams is not apparent. Middle: single photon point cloud for the 532 nm parallel channel on beam 3 
with the stand location indicated by the boxes. Bottom: beam 3 height distributions for the three stands showing 
number of received photons, binned at 0.5 m, with each channel normalized by its peak amplitude (left); foil 
cumulative distributions (FCDs) of total received photons, excluding the laser after-pulse, with FIOME indicated 
showing its height range where the channel distributions diverge (center); and canopy cumulative distributions (CCDs) 
of received photons excluding the ground return (right). Colors represent the 1064 nm perpendicular (orange) and 
parallel (red) channels, and the 532 nm perpendicular (light green) and parallel (dark green) channels. 






Fig. 4. Results for the Cedar Bridge Pine Barrens, NJ data acquired on 08/04/10. Top: SIMPL profile track (red) displayed 
on an aerial photograph from Google Barth acquired on 12/30/02. The stand segments are -100 m long. See Figure 3 
caption for explanation. 

clear cut. Seed trees left standing to re-populate the site are observed in the point cloud, attesting to the high horizontal 
resolution achieved by this small laser footprint, single photon data. 

In order to evaluate canopy vertical structure range histograms were computed by aggregating all detected photons in 
each of the 9 selected stands for each of the four channels on beam 3. The range histograms are presented in three ways 
in Figures 3-5, height distributions normalized by peak amplitude, full cumulative distributions (FCD) including the 
canopy and ground returns (excluding the laser after-pulse), and canopy cumulative distributions (CCD) (excluding the 
ground return). HOME is indicated on the FCDs, determined in the same manner as for conventional lidar waveforms. 
Where the distributions for the four channels at a site diverge, the height range of HOME is shown. The primary 
features of the individual stand height distributions for all four channels are very similar including: (1) the stand height, 
(2) the slope of the leading signal up to the initial canopy peak, where a greater slope corresponds to greater openness of 
the outer canopy surface 141 , (3) the crown depth indicated by the height range of the broad, primary canopy signal, (4) the 
location of prominent layers in the crown and the under-story indicated by narrow, secondary peaks, (5) the summed 
amplitude of the ground return relative to that of the total return which is the reflectance-weighted measure of canopy 
cover 141 , (6) segments in the under-story lacking significant cover (canopy vertical gaps) where the signal approaches 
zero, and (7) the width of the ground return which is a measure of the relief within the stand due to slope and/or 
roughness 171 . The consistency of the narrow peaks between all four channels for a stand indicates that these are real 



attributes of the stand structure, not noise, and demonstrates the high vertical resolution of the micropulse, photon 
counting technique even at this relatively coarse binning of 0.5 m. Differences between these height distribution 
attributes document significant structural differences between the stands. 

The similarity of a stand’s height distributions for the four channels leads to significant similarity in that stand’s FCDs, 
especially so for SERC stand A and HNF stand A and B, yielding a common HOME for all channels. The other stands 
show some divergence in the FCDs to varying degrees, with channels displaced to the right having a greater amount of 
received photons occurring above a given height. For SERC stand C and PB stands B and C the 1064 nm || channel is 
left-most and the 532 nm _|_ channel is displaced most to the right. For HNF stand C 532 nm || is left-most and 1064 nm 
_|_ is displaced most to the right. For the other stands the separations are less well defined. The divergence between 
channels leads to within-stand differences in HOME up to 6 m (SERC stand B) in an absolute sense and up to -20% as a 
fraction of stand height (SERC stand C and PB stand B). Although HNF stand C has a significant divergence of the 
channels, HOME is constant because the separation occurs at cumulative percentages greater than 50%. At some 
percentiles the corresponding heights can differ substantially. The largest case is HNF stand C where the 532 nm || 70 th 
% height is essentially zero and is 9 m for 1064 nm which is -50% of the stand height. 
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Fig. 5. Results for the Huron National Forest, MI data acquired on 07/29/10. Top: SIMPL profile track (red) displayed on a 
mosaic of nadir-viewing video frames acquired in-flight. The stand segments are -100 m long. See Figure 3 caption 
for explanation. 



On the FCDs the reflectance-weighted canopy cover is indicated by the percentile where the distributions reach the start 
of the ground return height. Divergence between the channels for a stand results in between-wavelength differences in 
the cover estimate that are for the most part small, below 10%. SERC stand C exhibits the largest difference, with 55% 
cover for the 1064 nm channels and 70% cover for the 532 nm J_ channel. Differences in the FCD intercept with the 
ground return document significant between-stand differences in their cover. For example, the three stands at PB do not 
differ markedly in appearance in the point cloud plot because their vertical structure is fairly similar, as indicated by the 
shape of their height distributions, but the cover ranges from 70% (stand A) to 45% (stand C). 

The generally close agreement between the 1064 nm and 532 nm HOME and canopy cover metrics documents that these 
canopy parameters, at least for these stands, are rather insensitive to wavelength differences, and indicates that biomass 
estimation using N1R and green laser altimeters will be comparable. This suggests that the 532 nm ICESat-2 
observations can be interpreted using similar biomass estimation relationships that have been developed for the 1064 nm 
SLICER, LVIS and GLAS instruments. 

In all the cases but one the within-stand channels CCDs are essentially the same. For SERC stand B there is some 
separation of the channels at high cumulative percentiles. Because the ground return signal is removed from the CCDs, 
they are a more direct means than the FCDs to compare canopy vertical structure, independent of cover. For example, 
the PB CCDs for stands A and B are very similar despite their differences in cover, consistent with their similar crown 
shapes in the height distributions. 

2.3 Discussion 

This work is the first of its kind, examining dual-wavelength, polarimetric lidar observations with high spatial and 
vertical resolution acquired through a diversity of canopies. As such, the results could not be fully anticipated. 
However, we had several expectations. We assumed that increasing multiple scattering with depth into canopies would 
cause an increasing amount of _|_ received energy as compared to ||, similar to water column volume scattering observed 
in SIMPL data 1231 . Therefore, we expected to see divergence between these channels in the CCDs progressing into the 
canopies with the || signal displaced to the right indicating that is preferentially occurring higher in the canopy than the 
_L signal. Furthermore, increasing reflectance generally increases the degree of multiple scattering so we expected that 
the effect would be more pronounced at 1064 nm due to the substantially higher foliage reflectance in the NIR compared 
to the green, often by as much as a factor of lOx. Therefore, the similarity across channels in the CCDs is unexpected. In 
the one case where divergence is seen (SERC stand B) the 532 nm result confonns to our expectation but the opposite 
relationship holds true for 1064 nm. Several factors could be diminishing detection of a depolarization signal with depth 
into the canopies. The laser footprint is very small, < 0.2 m. If it is small with respect to gaps in the canopy the first 
encounter with a plant surface, where single scattering dominates, could occur deep within the canopy thereby reducing 
the occurrence of multiple scattering at lower depths. In addition, the very small detector IFOV used to limit the solar 
background noise rate could be excluding a substantial portion of the multiply scattered photons. With decreasing IFOV 
the return signal becomes increasingly dominated by the single-scattered direct backscatter from the first encountered 
surfaces. The degree to which multiply scattered photons are excluded depends on the path length between scattering 
events, controlled by foliage shape and spacing. To better understand these unexpected results we are initiating 
modeling work to evaluate the effect of instrument parameters on the detection of canopy depolarization. 

Because the CCDs are insensitive to channel differences, we interpret the divergence of channels within the FCDs to be 
due to changes in the relative magnitudes of the canopy and ground retro-reflectance between channels. For the most 
common case of the 532nm _|_ channel being most displace to the right, that indicates that the ground is least reflective 
relative to the canopy for this wavelength and polarization state. The degree of channel separation for a stand is a 
measure of the magnitude of ground versus canopy reflectance differences. These differences in the reflectance ratios 
provides information on the optical scattering properties of the canopy and ground cover. The changes in the order of 
channel displacements between stands indicates there are changes in canopy and/or ground cover type. We will 
investigate the potential use of this information for differentiating canopy and ground cover types as a part of our 
modeling work. 

The number of stands we have assessed in this work is limited to nine. Further work examining more examples of these 
stand types at the sites is needed to ascertain if the relationships observed are consistent. In addition, the stands 
examined are mostly rather open, pine-dominated forest cover with only two closed-canopy, broadleaf forests included 
in the assessment. More broadleaf forest cover types with a greater diversity of stand structures will be included in 
future work. Lastly, the wavelength comparison of HOME and cover was done using micropulse, photon counting for 



each. We do not expect that results will differ substantially between conventional waveforms and photon range 
histograms because both are observing the travel time of photons in an unbiased manner. However, this assumption 
needs to be confirmed by direct comparison of NIR waveforms and green photon counting height distributions. 


3. CONCLUSION 

This work provides the first detailed examination of forest canopy structure as a function of two laser altimeter 
wavelengths and two polarization states. Analysis of results for nine stands with different age, structure and tree 
composition at three sites in Maryland, New Jersey and Michigan shows that the SIMPL instrument’s implementation of 
the micropulse, photon counting approach provides structure measurements with high vertical and spatial resolution. 
The well-resolved forest cover structure portrayed in single photon point clouds enables selection of distinct stands for 
subsequent analysis using height distributions of received photons analogous to conventional lidar waveforms. The 
primary attributes of the examined stands are similar at 1064 nm and 532 nm, in particular for the HOME and cover 
metrics commonly used for above-ground biomass estimation, indicating that NIR and green lidar systems will provide 
comparable estimations. Significant differences expected in height distributions as a function of wavelength and 
polarization state, due to within-canopy volume scattering causing increasing depolarization with depth, were observed 
for only one stand. We speculate that the small lidar footprint and/or small detector IFOV are minimizing detection of 
multiply-scattered photons with polarization perpendicular to that of the transmitted laser pulses. Modeling work is being 
initiated to assess instrument parameter influence on the detection of multiple-scattering. We infer differences in canopy 
cover estimates as a function of wavelength and polarization are related to differences in the ground reflectance relative 
to that of the canopy between the four channels. Our modeling work will assess the potential for differentiation of 
canopy and ground types based on the differences in these reflectance ratios. This work is providing pathfinder results to 
better understand and interpret the non-polarimetric, 532 nm micropulse, photon counting data that will be acquired by 
the NASA ICESat-2 mission, scheduled for launch in 2016. 
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